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Using a continuum Dirac theory, we study the density and spin response of zigzag edge terminated 
graphene ribbons subjected to edge potentials and Zeeman fields. Our analytical calculations of the 
density and spin responses of the closed system (fixed particle number) to the static edge fields, 
show a highly nonlinear Weber-Fechner type behavior where the response depends logarithmically 
on the edge potential. The dependence of the response on the size of the system (e.g. width of a 
nanoribbon) is also uncovered. Zigzag edge graphene nanoribbons, therefore, provide a realization of 
response of organs such as the eye and ear that obey Weber-Fechner law. We validate our analytical 
results with tight binding calculations. These results are crucial in understanding important effects 
of electron-electron interactions in graphene nanoribbons such as edge magnetism etc., and also 
suggest possibilities for device applications of graphene nanoribbons. 
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I. INTRODUCTION 

Since it's first isolation^, graphene, a single layer of 
carbon atoms arranged on a honeycomb lattice, has 
attracted much attention. Several unconventional and 
novel phenomena have been observed in graphene; these 
include quantum Hall effect at room temperature, Klein 
paradox etc., to name a few^— Many of its interesting 
properties originate from the fact that energy spectrum 
of graphene, a zero band-gap semiconductor with lin- 
early vanishing density of states(DOS) near the chem- 
ical potential, resembles the Dirac spectrum of mass- 
less Fermions£ Due to it's unusual electronic proper- 
ties, graphene holds promise of replacing semiconductor 
based transistors in future nanoelectronic and spintronic 
devices £ Advent of several electronic devices, such as 
single electron transistor based on graphene flakes^ and 
nanodots of size as small as ~ 30 nmp quantum interfer- 
ence devices^i etc. are just the beginning of a fascinat- 
ing future. 

Several interesting nanostructures can be fabricated 
from graphene. The simplest of them is the quasi 
one-dimensional nanoribbon, made by terminating the 
graphene sheet by two edges of either "armchair" or 
"zigzag" type^i Such terminations can significantly in- 
fluence the electronic structure of the graphene nanorib- 
bon. In particular, tight binding analysis of graphene 
bound by two zigzag edges shows that there are local- 
ized electronic states at the edge with nearly flat en- 
ergy dispersion^! In terms of magnetic structure, half 
filled zigzag edge nanoribbons are reported^ to have an- 
tiferromagnetic ground state, with very large moments 
localized at the edge sites compared to bulk and the 
magnetizations of the opposite edges being antifcrromag- 
netically aligned. The investigation has further been 
extended to zero dimensional graphene nanodots (half 
filled or undoped) of different shapes such as rectan- 



gular, hexagonal etc . 13 ' 14 , which are also reported to 
be antifcrromagnctic (zero net magnetization) above a 
critical size. There has also been interest in valley 
physic o 15 ! 16 , giant magnetoresistanc o 1 7 ' 18 phenomena in 
graphene nanoribbons. These developments strongly 
motivate studies aimed at understanding the electronic 
properties of graphene nanosystems and study of their 
response to stimulus. 

Responses of quantum systems, including nanosys- 
tems, arc usually linear in the stimulus, with response 
functions described by Kubo formulae^ Nature also 
provides realization of systems that inherently have a 
strongly nonlinear response. Key examples of these are 
the eyes and ears of living organisms whose response is 
described by the Weber-Fechner lawi2£ If S is the stimu- 
lus and R is the response, a system obeying the Weber- 
Fechner law satisfies 

where AR is the change in the response for a change AS 
of the stimulus. This implies that the response depends 
logarithmically on the stimulus. Indeed, this is why eyes 
and ears of living organisms can sense light and sound of 
intensities that range over 10-15 decades! Indeed, it is 
interesting to explore quantum/nano-systems that have 
such a nonlinear response. 

In this paper we demonstrate theoretically that zigzag 
edge terminated graphene nanoribbons subjected to edge 
potentials show responses that are highly nonlinear and, 
most interestingly, obey the Weber-Fechner law. Our 
motivation for this work arose from the earlier works 
on edge state magnetism cited above that we found to 
be very robust 31 Edge potentials can therefore be ap- 
plied externally or can be "self generated" by the system 
due to interactions (as in the case of edge state mag- 
netism), and the considerations of this paper apply to 
both cases. Specifically, we analytically solve, using a 
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continuum Dirac theory for graphcne, for responses such 
as particle density and spin density (of the closed system 
with fixed number of particles) to applied static edge 
potentials and Zeeman fields in zigzag edge terminated 
graphcne nanoribbons. We show that they have a Weber- 
Fechner type nonlinear response that can be tuned using 
the size (width) of the ribbon (see eqns. (|5^|) and (|53[l). 
Indeed this work provides an uncommon example of a 
quantum system whose nonlinear response can be calcu- 
lated analytically, in addition to suggesting many novel 
uses of zigzag terminated graphene nanoribbons. We also 
compare our analytical calculations with full tight bind- 
ing simulations and demonstrate excellent agreement be- 
tween the two. 

The paper is organized as follows. Section |TT] reviews 
the continuum field theory of graphene and the earlier 
work on the electronic structure of zigzag terminated 
nanoribbons, and in addition introduces the notation 
used in the paper. In Section IIIII we state and solve 
the problem of the zigzag terminated graphene nanorib- 
bons subjected to edge potentials. This is followed by 
sections IIVI and [V] where the density and spin responses 
to the edge potential are discussed and shown to have 
Weber-Fechner type behavior. These analytical results 
are shown to be in excellent agreement with full tight 
binding calculations in Section IVI1 The paper is con- 
cluded in Section IVHI which includes a discussion of pos- 
sible applications and future directions of research. 
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FIG. 1. (color online) Graphene structure/zigzag edge termi- 
nated nanoribbon-geometry explaining notation used in the 
text. The triangular lattice parameter and the "unit cell 
width" of the nanoribbon is denoted by a; the atoms enclosed 
by the two dashed lines form the repeat unit basis along the 
^-direction. W is the width of the nanoribbon. The two sub- 
lattices A and B are color coded in red and blue respectively. 
The hopping amplitude for electrons is denoted by t. 



II. GRAPHENE NANORIBBONS 

The purpose of this section is two fold: to present a 
brief review of the electronic structure of graphene and 
graphene nanoribbons (for a detailed review cfi^), and in 
the process introduce the notation used in the text. 

The low energy electronic degrees of freedom in 
graphene are the p z orbitals on each carbon atom. Elec- 
tron hopping, with amplitude t, from p z orbital of one 
carbon atom to one of its three neighbors in a honey- 
comb lattice (sec Fig. [T]) is described by a tight binding 
model 

Utb = -tJ2 c l c i° (2) 

where c\ a (ci a ) creates (destroys) an electron of spin a at 
lattice site i. As is well knowni, diagonalization of this 
hamiltonian leads to two bands which "touch" at two 
non-equivalent points K + and K in the first Brillouin 
zone of the honeycomb lattice - called the Dirac points - 
that are related by time-reversal symmetry. By particle- 
hole symmetry of the hamiltonian, the energy where the 
bands touch is zero. 

Each carbon atom contributes one p z electron, so the 
chemical potential in undoped graphene is zero and the 
highest occupied states are at the Dirac points. For 
this case, and for lightly doped graphene, the electronic 



excitations can be conveniently described by a contin- 
uum theory. The electronic low energy excitations (for 
each spin a) are described by two component wavefunc- 
tion ip aa (r), where a denotes the "sublattice flavor" 
and can be A or B (see Fig. [TJ. Moreover, since the 
Bloch wavefunctions that contribute to these wavefunc- 
tions mainly arise from the neighborhood of K + and K 
points (called as the + and — valleys), the wavefunction 
can be expanded as 

Vw(r) = ^ aa+ (r)e iK +- r +^ a .{r)e iK -- r (3) 

where ip aav (r) (v = ±), stands for wavefunctions with 
"valley flavor". For each spin, therefore, the wavefunc- 
tion is described by a four component spinor, ^ a (r) = 

(-0 CT A+( r ) VVB+(r) ipaA-(r) ipaB-{r)) T , which satis- 
fies the Schrodinger equation 

U K ^ a = EV a (4) 
with E the energy eigenvalue, where 

Uk = hv F (^ (T Q P _ J? p J (5) 

where vp = ^ is the Fermi velocity at the Dirac points, 
P = — t/iV and cr = e x o x + e y a y ; a x y are the Pauli ma- 
trices. The energy eo is related to the hopping parameter 
via eo = : ^ L - The hamiltonian has an emergent SU(4) 
symmetry^ The eigenstates of the hamiltonian are 
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FIG. 2. (color online) Brillouin zone of graphene. The points 
K± = (±§^,0) are the Dirac points. The line segment 
Q Q+ represents the Brillouin zone of a zigzag nanoribbon 
(NRBZ) (segments in blue and red) with Q± = (±f,0). 
Segments E+Q+ and Q E corresponds to "edge Brillouin 
zone" (EBZ) (segments in blue) and is one third of the NRBZ 
(E± = (±|^,0)). For the nanoribbon, the Dirac point Jf_ is 
equivalent to E+, and K+ to E-. 



linearly dispersing modes with energy ?wf|A;| about each 
of the Dirac valleys. 

Consider, now, an external potential V a (r) that acts 
on the lattice. Assuming that the potential is "slowly 
varying" , the Hamiltonian becomes 



H=H K +Uv 



where 



H v = 



V 
V 



V = 



V A (r) 
V B (r) 



(0) 



(7) 



and states evolve via the Schrodinger equation (|4]) with 
V. instead of Hk: ■ 



A. Review of earlier work: key results 

We now briefly review the results of Brey and 
Freti g 23 ' 24 for states in a graphene nanoribbon that are 
terminated by zig-zag edges in the absence of any ex- 
ternal potentials. Edge states including spin-orbit inter- 
action have also been studied^ but we shall not con- 
sider this here. The nanoribbon occupies the region 
— oo < x < co, < y < W as shown in Fig. [TJ with 
the edge at y = W being terminated by A type atoms 
and that at y = by B type atoms. The nanoribbon 
may also be thought of as a one dimensional crystal with 
a basis of atoms as shown in Fig.[T] The Id Brillouin zone 
of this system, shown in Fig. [2] is the segment Q Q + . 
Within this scheme, the Dirac point K_ is equivalent to 
the point E + in the Id Brillouin zone, and similarly K + 
is equivalent to 
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FIG. 3. Energy gap of the edge-states as a function of the 
momentum in EBZ. 



The problem has translational symmetry in the in- 
direction, and the energy eigcnstatcs of §5§ are of the 
form i£ a (x,y) = ip aa v(x,y) = e lkx 4> aa v{y), where k is 
a wave-vector (momentum) measured from the Dirac 



points (E+, EJ). The boundary conditions^ are 
VfeA±(0) = 0; i>„ B ±{W) = 0. 



(8) 



The equations for the valleys therefore nicely decouple, 
i.e., for zigzag edge nanoribbons, valley index is a good 
quantum number. 

We focus here on the edge-states which play a piv- 
otal part in our work. Edge-states are those eigenstates 
for which 4> aa v(y) are appropriately "exponentially de- 
caying" functions, i.e., they are wave propagating along 
the a;-direction with amplitudes "localized" on the edge 
atoms. Thus, for a given k, 



4>aav(y) = F„ 



G 



(9) 



q = sjk 2 — (e/aeo) 2 where e is the energy eigenvalue, 
i.e., %k.^(j — £^a, F and G are constants. One finds 
that q satisfies the equation 



tanh (qW) 



(k 



(10) 



where ( = —v is +1 for E+ (K- valley) and —1 for E- 
(K+ valley). It is easily concluded that edge-states are 
possible only for momentum in the "edge Brillouin zone" 
(EBZ, see Fig. [2J [Q_, E_]\J[E + , Q + ], and there are only 
"right movers" at the E + valley and, only "left movers" 
at the E_ valley. In particular, there are edge-states 
around E + only for k g [k c , j^], and for k 6 [— — k c ] 
around E_, where 



1 

W 



(11) 



Note that in the limit of W ^> a, we see that EBZ is about 
one third of the BZ of the nanoribbon, i.e., roughly, there 
is one edge-state for every three edge atoms. Now, for 
each k in the EBZ, there are two edge-states with energy 



4 



±e(fc). For large k (~ 7r/3a), one obtains the energy 
dispersion as 



eo 



±2kae 



-kW 



£g( fc ) 
eo 



4fcae 



-kW 



(12) 



where the + sign is for the "conduction band" , and the 
— is for the valance band. One can define a gap function 
as the difference in the energies of the conduction and 
valance band for each k. It is interesting to note the re- 
markable exponential distribution of the energy gap (see, 
for example, Fig. [3J); as will become evident, this charac- 
ter of zigzag edge terminated nanoribbons play a crucial 
role in their nonlinear response. The edge-state wave- 
functions, to be used in the discussion below, are 



<P<ta-(v) 



■2<1 



sinh (2qW) - 2qW 



smh(qy) 



(13) 



tB-{v) = T\ 



Bwh(2qW) - 2qW 



smh(q(W — y)) 



where the =F stands for valance and conduction bands. 
These wavcfunctions arc for the E+ (K-) valley, with 
similar expressions for those at E_ . These analytical re- 
sults have been tested against full numerical calculations 
of the tight binding model.— 

In addition to the edge-states, the nanoribbon also sup- 
ports "bulk states" and these can be obtained by letting 
q in (fit)]) to be imaginary. The key point to be noted 
is that in a finite width nanoribbon, the energy of the 
bulk states, in magnitude, are at least ~ aeo/W - in this 
sense the bulk states are "fully gapped" for narrow rib- 
bons. The same argument applies to energy differences 
between different bulk bands, they are of the order of 
- aeo/W. 



III. GRAPHENE NANORIBBONS WITH AN 
EDGE POTENTIAL 



The specific problem of interest to us is when the 
nanoribbon is subjected to edge potentials. In terms of 
the potential operator V, the edge potential ([7]) can be 
expressed as 



V 



aV A 5(y -W) 

aV B S(y) 



(14) 



where S(-) is the Dirac delta function, Va and Vb have 
units of energy, and represent, respectively, the values of 
constant potentials acting on the edge atoms A (y = W) 
and B (y = 0) (see Fig. [I). It is also possible to generalize 
the problem to include edge Zceman fields as will be done 
later in the paper. 



A. Qualitative Discussion 

We first discuss why such an edge potential flT4l pro- 
duces highly non-linear Weber-Fechner like responses. 
Note that the potential (fbfj) retains the translational 
symmetry of the problem in the x-direction, i. e., the mo- 
mentum k introduced in the previous section continues 
to be a good quantum number. Also, the potential does 
not mix valley states. Therefore, the eigenstates on the 
application of the edge potential will be linear combina- 
tion of states with the same momentum k. Consider now 
\Va\ ~ \Vb\ ~ V, where V is a characteristic magnitude 
of the edge potential with V <C aeo/W < eo- Based on 
the discussion in the previous section, it is evident that 
the applied edge potential will have little effect on the 
bulk states due to large energy denominators that result 
in the quantum mechanical mixing by the edge potential. 
Edge-states are however not spared. Some of them are 
very strongly affected by the edge potential. In particu- 
lar states with k > k m are strongly affected in that the 
edge potential mixes the valance and conduction band 
edge-states. We can estimate k m as follows: the matrix 
element of V that mixes the conduction band \k^// 3e ) and 



valance band states \k^f ge ) is roughly (using (fT3")) ) 



j(k) = m 



edge 



V\k 



edge\ 



AkaV. 



(15) 



Strong mixing of states is obtained whenever 



e g {k)<l{k) 



Akae e- kw < 4kaV. (16) 



Thus we find that all the edge-states with k > k m are 
strongly affected by the edge potential, with 



k m = -T77I11 ( — 
W \e 



(17) 



and those with k < k m are relatively unaffected. In- 
terestingly, we do not expect any response when V is 
below a threshold V t h = eoe - "^. The physical origin 
of the Weber-Fechner like responses is now immediate. 
The response of a system is typically a change in the 
charge or spin density from the case with no edge po- 
tential. It is clear that the mixing of conduction and 
valance band states can push the valance band edge- 
states with k > k m above (or bring conduction band 
edge-states below) the chemical potential. The response 
to the edge potential arises solely from the change in oc- 
cupation of these states. Other edge-states k < k m and 
bulk states are essentially unaffected and do not con- 
tribute to the response. These simple arguments sug- 
gest that such edge potentials will evoke highly nonlin- 
ear response of the Weber-Fechner kind in a zigzag edge 
terminated graphene nanoribbon. Observe also that k m 
depends inversely on the width W of the nanoribbon. 
Motivated by this discussion, we shall now obtain the 
edge-state dispersion and wavefunctions in presence of 
the edge potential of the form flT4"|) to make these ideas 
precise. 
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B. Edge-states in presence of an edge potential 



As noted, the edge potential (fT4f preserves the transla- 
tional symmetry in the x-direction and the energy eigen- 
states have the form f/Va- y) 
E+(K-) valley, and satisfy 



0.010 



0.005 



' lkx 4>au-(y) for the 
e a{D y + k)4> aB -(y) + V A aS(y - W)4> aA -{y) = £<f>crA- 



(y) or 

(18) 



.0.000 



e a(D y - k)<j>o-A-(y) + V B a5{y)(j> a B-{y) = s<j> aB -{y) 

(19) 



-0.005 



where e is the energy eigenvalue and D y is the derivative 
with respect to y. The equations (jT5J) and (jT9")l can be 
solved in the region < y < W as 



b aA .(y) = Fe™ + Ge 
1 



(20) 



^b- (y) = - (F(q - fc)e w - G(q + k)e~™) (21) 

£ 

Integrating (fT5| from W — r\ to W + r\ and letting r\ — s- 0, 
we obtain 

e ^B- (W0 + Va0*a- (W) = (22) 
and by integrating ()19j) from y = —77 to y = r/, we get 

eo^rA-(O) +Vb&,b- (0) = (23) 



Equations (|22j) and (|23|) provide the necessary boundary 
conditions. The resulting secular equation for the energy 
eigenvalue reads 



'i 



e% J tanh(gVF) 



('(> 



e 
ae 



VaVb \ 



1 ^~) k (24) 



The secular equation for the E_(K + ) valley will have a 
— sign on the right hand side. 



C. Half Sheet 

We shall first consider a half sheet of graphenc which 
occupies the region y > - the edge is terminated by 
B type atoms (see Fig. [1} and a potential V r (<C eo) acts 
on these edge atoms. The secular equation (|24p can be 
readily solved to obtain the dispersion 



e(fc) 



2eoV_ 



ak 



V 

2—ak, 

eo 



(25) 



near E + (a negative sign is obtained for states near E_ ) . 
A constant edge potential on a half sheet makes the edge 
modes linearly dispersing, i. e., moving with a constant 
velocity which is proportional to the applied potential. 
As anticipated in the qualitative discussion all the edge 
states are affected (since k m wO). Moreover, these edge- 
states are such that, for a positive edge potential V, the 
states near E + are right movers, and those at E_ are left 
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FIG. 4. (color online) Dispersion of edge states in the presence 
of a symmetric edge potential. Blue line is the valance band 
and red line the conduction band. Solid lines are for a positive 
value of the edge potential, while the dashed ones are for a 
negative value of the edge potential. Note the near linear 
dispersion of the edge states near Q. 



movers. Most interestingly, by applying a negative edge 
potential, the sign of the velocity can be reversed, i. e., 
states near E + become left movers etc.! The edge-state 
wavefunctions (at E+) are 



<f>*B-(y) 



2q 



(1 + (V/e ) 



!>*A-(y) = - V(j> aB -{y) 

(26) 



where q 

It must be noted that for the half sheet, the bulk states 
near E + are also affected by the edge potential. However, 
the bulk dispersion (and gap) is not exponential as in the 
case of edge states. Therefore, only small corrections to 
the bulk states are expected in a small region near the 
Dirac points. 

We now proceed to consider the case of finite width 
nanoribbons. We shall consider two special cases of edge 
potentials. First is a symmetric edge potential where 
Va = Vb = V and second is an anti-symmetric edge 
potential V A = —Vb = V. 



D. Nanoribbons with Symmetric Edge Potential 

A symmetric edge potential V A = V B = V, shifts the 
energies of the conduction and valance bands in the same 
direction, i. e, for every k, the energies of both the bands 
cither increase or decrease depending on the sign of V. 
For such a potential with V A = V B = V, the secular 
equation (|24p reduces to 



eg J tanh (qW) 



! I — 



= 1 



ae 



V 2 \ 
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where e, the energy eigenvalue, is one of ±eo<J\/fc 2 — q 2 . 
The two solutions evolve from the conduction (valance) 
bands of the zero edge potential problem. For the case of 
the valance band, we see that edge state solutions appear 
only when k > k^ , where 



' 



V 1 



(eo + V) 2 W 



(28) 



and similarly for the conduction band only when k > fcj? 
with 



V 2 1 



(e - V) 2 W 



(29) 



For edge potentials with V -C eo, k^ and k^ are only very 
slightly different from the case without the edge potential 
and do not affect the physics in any significant way. 
When k is near k c , the dispersion is essentially unaffected 
by the presence of a small edge potential. However, for 
ks close to the EBZ boundary, k sa the dispersion for 
both the conduction and valance bands is 



:(k) 2V 



eo 



ak 



eo 



(30) 



much like the case of the half sheet discussed in the pre- 
vious section. Fig. @] shows both the bands in the EBZ 
obtained from the numerical solution of ([27)) ; they have 
all the features anticipated in the qualitative discussion 
above - states near the EBZ center (near E) are un- 
affected while those near EBZ boundary (near Q) are 
strongly affected. An important feature is that for a 
positive (negative) edge potential the valance (conduc- 
tion) band crosses the zero energy, and the conduction 
(valance) band has a minimum (maximum). It is clear 
that the wave vector k s m (the superscript s indicates that 
we are dealing with a symmetric edge potential) at which 
the valance (conduction) band crosses the zero of energy 
has a key role to play in the response. This wavevector 
k^ can be immediately obtained from (|2T|) : 



1 

W 



In 



(31) 



Near k m , the band that crosses the zero of energy dis 
perses linearly with a velocity v m given by 

'\V\ 
eo 



Vm 
Vf 



2e V 



-V 2 



hi 



(32) 



which we sec is independent of the ribbon width. 

Due to the symmetry of the edge potential the wave- 
functions are such that the total weight on the A sublat- 
tice is equal to the total weight on the B sublattice. 



E. Nanoribbons with Anti-Symmetric Edge 
Potential 



In the second case, we consider an antisymmetric edge 
potential such that Va = —Vb = V. This edge potential 
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FIG. 5. (color online) Dispersion of edge states in the presence 
of an anti-symmetric edge potential. Blue line is the valance 
band and red line the conduction band. Note the near linear 
dispersion of the edge states near Q; also, the lowest energy 
gap obtains at kg. 



is such that it mixes the conduction and valance band 
states, with a resulting "level repulsion". The secular 
equation (f2"4")l now becomes 



1 



V- 



e\ ) tanh (qW) 



1 



v- 



(33) 



for states near E+ (a negative sign on the r. h. s. will 
give the secular equation for states at E_), Moreover, 
in this case, energy eigenvalues come in ± pairs due to 
the (anti-)symmetry of the potential. We find that edge- 
states appear whenever fc > fc", 



k c . — 



V 2 1 



V 2 w 



(34) 



Again, we find that that states with k near k c are un- 
affected by a small edge potential, while those near the 
EBZ boundary arc strongly affected and disperse linearly 



£(fc) 

eo 



2V 
=F — ak 

eo 



(35) 



for the valance and conduction bands respectively. Fig. [3] 
shows the dispersion of the valance and conduction bands 
obtained from the solution of (|33]l . The smallest gap is 
no longer at the boundary of the EBZ, but occurs at a 



wavevector kg which is calculated as 



k a 
K g 



W 



In 



rn 

eo 



In - In 



\V\ 
eo 



and the energy gap is 



eo 



eo 



W e 



(36) 



(37) 



with, again, a strongly non linear dependence on the ap- 
plied edge potential. The wavefunctions of this problem 
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have qualitatively different features from the case of the 
symmetric potential. The wavefunctions are relatively 
unaffected for states with k sa in that the weight in 
the A sublattice and B sublattice are about equal. How- 
ever, assuming V > 0, most of the weight is in the B 
(A) sublattice for valance (conduction) band states with 
k > k^ . In fact the difference of weights is a monotonic 
function of k (see Fig. [BJ to be discussed in detail later). 

Armed with the discussion of the eigenstates of the 
system in presence of edge potentials we now proceed to 
discuss the responses. Results similar to those in fig. [4] 
and fig. [5] have been obtained by Wimmer— . 

IV. DENSITY RESPONSE WITH EDGE 
POTENTIALS 

This section focusscs on the zero temperature den- 
sity response of zigzag edge terminated graphene systems 
that are subjected to an edge potential. Throughout this 
section, we shall assume that the chemical potential of 
the system is maintained at zero. We shall discuss sepa- 
rately the three cases that were discussed in the previous 
section. 



A. Half Sheet 

From the dispersion for the edge-states (|2~5l) . it is im- 
mediate that for a negative edge potential, the edge states 
are pushed below the zero chemical potential. For a small 
edge potential, there is no further change of occupancy, 
and therefore there is no density response. However, for a 
positive edge potential, the edge states are thrown above 
the chemical potential, and therefore every edge-state 
in the EBZ is left unoccupied. It is evident, therefore, 
change in number of electrons (per repeat unit in the x 
direction) is 

AN = AN t + AN, = -^B(V) (38) 

o 

where 0(-) is the Hcaviside step function, the factor 3 in 
the denominator arises from the fact that the EBZ is a 
third of the full ID BZ, and the factor 2 in the numerator 
from two spins. It is also interesting to study how this 
loss of particles is distributed on the half sheet: 

am*) = -WrAff A (i - e-^d + ^)) 

(39) 
(40) 

fe 2 — V 2 ) 

where v = ^ e Sj + y 2 j ■ It is interesting to note that ha (y) + 
%(y) has an integrable power law decay, with diverging 
higher moments. 




FIG. 6. (color online) Contribution to the "charge density 
wave" amplitude from edge modes. The total "charge den- 
sity wave" amplitude C is well approximated by p apP rox{k) = 
0(fc — km) shown by the blue line; the red line is the exact 
result. 

B. Nanoribbons with Symmetric Edge Potential 

The change in the number of particles in the case can 
be calculated as follows. Note from Fig. |4] that the occu- 
pancy of all states with k < k s m is unaffected by the edge 
potential. For the case with V > (V < 0) and k > fcf„, 
the valance (conduction) band states are lifted above 
(pushed below)the chemical potential and are hence un- 
occupied (occupied). One finds that the change in the 
number of particles per repeat unit is 

AN = AN, + AN, = 2 sgn (V) Q + JL m (M) ) 

(41) 

which has a strongly nonlinear Weber-Fechner like be- 
haviour. It must also be noted that, the response will 
vanish for very small edge potentials. The threshold po- 
tential for a nonzero response is obtained as 

7T W 

V th = e e-—. (42) 
We can rewrite (|4T|) as 

AA/.-2sgn(y)-^ln(M), (43) 

which is in the standard Weber-Fechner form. 



C. Nanoribbons with Anti-symmetric Potential 

For an anti-symmetric edge potential, it is immediate 
that net particle number response vanishes for any V 
since the occupation numbers of the bands are unaffected. 
However, as noted in section PlI El the electronic states 
are such that the weight on the A sublattice is different 
from the B sublattice resulting in a "density wave" like 
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response. The density wave response associated with an 
edge-state of momentum k is 

p{k)= dy{\^ B -(y)\ 2 -\^ A -(y)\ 2 ) (44) 
Jo 

and this is a monotonic function of k. Fig. [5] shows a plot 
of p(k) as a function of k for a particular 1/ > 0. Now 
the "density wave order parameter" D is 

2a r^ a 

D = — dkp(k). (45) 

T Jo 

A detailed analysis shows that -D can be very accurately 
approximated by taking 



p(fc)">Pa(ft)=Q(*-C) 



where 



resulting in 



1 , (V3\V\\ 



(46) 
(47) 

(48) 



whenever \V\ is greater than an appropriate threshold 
as in (|42[) . For a positive edge potential, the "density 
wave" that arises will be such that density is higher on 
the edge B, than on edge A, i. e., a build up of parti- 
cles at the B edge at the expense of the A edge. This 
observation will be important in understanding the spin 
response of zigzag edge terminated graphene edge sys- 
tems as discussed in the next section. 



V. SPIN RESPONSE WITH EDGE ZEEMAN 
FIELDS 

In this section we will explore spin response of zigzag 
terminated graphene structures to applied edge Zeeman 
potentials. The edge Zeeman potential operator is 



V H * = 



-aaH A 6{y - W) 






-aaH B 5{y) 



(49) 



where Ha and Hb are Zeeman energies due to appropri- 
ate fields on edges A and B respectively. As noted in the 
introduction, such edge Zeeman potentials are encoun- 
tered in various mean field theories of graphene systems 
which treat interactions between electrons. We shall now 
use the results of the previous sections to obtain results 
for spin response in various cases. 



A. Half Sheet 

For the case of the half sheet with a Zeeman potential 
Hb = H, we find the dispersion (|25[) as 



e a {k) -2ae a H aH 

= —r —aft w — z ak 

e eg + H z e 



The spin response can be obtained by considering the 
filled states, and we obtain the magnetization per repeat 
unit in the x direction as 



M = AiV t - ANi = 
= I sgn (if) 



-l @ (-H)-(-±e(H)) 



(51) 



Not unexpectedly, even an infinitesimal Zeeman field pro- 
duces a large magnetic response with a magnetization of 
1/3. The spin density, again, dies as a power law simi- 
lar to (f3T)|). with a larger magnitude on the B sublatticc 
atoms and a smaller one of opposite sign on the A sub- 
latticc. 



B. Nanoribbon with Symmetric Zeeman Potential 



For this case Ha = Hb = H, one can show by a par- 
ticle hole transformation c] ja — > (—) a c Sa (a is ±1 for 
a = A, B) , the chemical potential is zero when the to- 
tal number of electrons is maintained at one electron per 
site. We have a situation where the total number elec- 
trons do not change, but the number of electrons of each 
spin species changes to produce a magnetic response. We 
find that the magnetization of the system is (see (|41|Q 



M = AN* - AN = 2 sgn (H ) ( - + — In ( S 



2) 



ak. (50) 



for a Zeeman field larger than a threshold H t h given in 
(|4"2")k again a Weber- Fechner like behaviour. Since the 
wavefunction of the edge-states have a significant amount 
of weight at the edges, we find that such a Zeeman edge- 
potential produces a magnetization that is primarily at 
the edges. In the present case the magnetizations of 
both the edges are parallel - a "ferro" configuration - 
with magnitude roughly given by half of (|5"2"j) . The band 
structure that emerges in the presence of the symmetric 
potential is such that, at the chemical potential, the dis- 
persion is linear for both spins. However they disperse 
with opposite velocities given by (|32p with V replaced 
by H. The direction of the velocities are reversed at the 
other valley. 



Nanoribbon with Antisymmetric Zeeman 
Potential 



For the antisymmetric Zeeman potential, the chemical 
potential again is zero (one electron per site). Since the 
total field is zero, there is neither a total density or to- 
tal spin response. Rather, such a Zeeman field produces 
a "spin density wave" like response where the magneti- 
zation on the A sublatticc is opposite of that on the B 
sublattice. For a Zeeman field larger than an appropri- 
ate threshold (see (|42p) the total magnetization on each 



9 



sublattice (per unit cell) is 



bons when interactions are included^. 



M S = -M, = Sg n(F)(i + ^ln(^))(53) 



D. Doped Nanoribbons with Edge Zeeman Fields 



with a Weber-Fechner like dependence on the edge mag- 
netic field (Zeeman energy). The resulting configuration 
is such that most of the spin polarization of each sublat- 
tice resides at the edges resulting in an edge magnetized 
state with an "anti-ferro" orientation of the edge mo- 
ments. The system is is gapped with a gap given by (|57|) 
with V replaced by H . 

It is particularly interesting to note that for a given 
magnitude H of the edge Zeeman potential, the magni- 
tude of the edge moment for the anti-symmetric potential 
is larger than that for the symmetric potential. This is 
key to understanding the magnetic properties of nanorib- 

I 



It is interesting to investigate the properties of doped 
nanoribbons, particularly with a view towards possibility 
of application in electronic devices. 



1. Symmetric Zeeman Potential 

We consider a zigzag edge terminated graphene 
nanoribbon with a small doping 8 per edge atom (<5 is 
a dimcnsionless number). A positive 5 corresponds to 
hole doping, and negative 6 to electron doping. 



As noted earlier, the dispersion near the zero energy (i.e., near k ss k m ) is given by 



cr(fc) / v m . . 

= -o — a(k - k m ) - 

e \vf 



W 

I s m a (k — k m ) 
a ' 



(54) 



where v m is as defined in (|32p with s r 



epH 



H 2 



4-h 2 



In 



H\ 



('a 




(55) 



is a dimcnsionless number that depends on the poten- 
tial. Using this dispersion relation, we can calculate the 
change in the magnetization per edge atom as 



AM{6) = - 




to lowest order in 5. We see that the edge magnetization 
falls quadratically with doping with coefficient propor- 
tional to the width of the ribbon. The change in energy 
(per repeat unit in the ^-direction) as a function of dop- 
ing is 

2. Antisymmetric Zeeman Potential 

In this case the magnetization per edge atom falls as 

AM = -\S\ (58) 

when the doping per edge atom is 5. Energy changes 
according to 

eo e 
where e g is the energy gap given in (|37l) . 



VI. COMPARISON WITH FULL TIGHT 
BINDING SIMULATIONS 

In this section we present comparison and validation of 
the theoretical results of the previous sections with full 
tight binding calculations of undoped and doped systems. 
We discuss only the case with the edge Zeeman fields. 
The tight binding calculations are performed using the 
momentum basis to represent @ with appropriately ap- 
plied edge potentials. We have used up to 10000 fc-points 
in the BZ Q Q + (see Fig. [5]) in the calculations reported 
here. 

Fig. [7] shows a plot of the calculated moment M per 
edge using the tight binding calculations as a function of 
the width W of the nanoribbon with a symmetric edge 
Zeeman field H. Results are shown for five orders of 
magnitude for the Zeeman field H. We find that the 
Weber-Fechner behaviour predicted in (|52j) is well repro- 
duced in the tight-binding calculations. In particular, 
the coefficient of the (a/W) term in (|5^|) is reproduced 
to with in a percent. 

Results of a similar study of nanoribbons with an an- 
tisymmetric Zeeman potential is shown in Fig. [51 where 
again we see that Weber-Fechner law predicted ([53)1 is 
realized. For a given H and W, we see that the anti- 
symmetric potential produces a larger edge moment, as 
predicted in and (|53"|) . 

We have also investigated various other quantities 
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FIG. 7. Dependence of edge magnetization M on the width W 
of the undoped graphene nanoribbon for different symmetric 
edge Zeeman fields. Points are results of the tight binding 
calculations, lines are fits to theory. 
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FIG. 8. Dependence of edge magnetization M on the width 
W of the undoped graphene nanoribbon for different anti- 
symmetric edge Zeeman fields. Points are results of the tight 
binding calculations, lines are fits to theory. 



(e.g., the gap (|5T|)) using the tight binding calculations 
and compared them with our analytical predictions and 
found excellent agreement in all cases 2 ^. 

The dependence of the magnetization and energy on 
the doping 8 obtained from tight binding calculations is 
shown in fig.O The results (JS7J), d£5J} and O are 
recovered from the tight binding calculations both for 
hole and electron doped systems. 



VII. SUMMARY AND CONCLUSION 

In summary, we have shown that zigzag edge ter- 
minated graphene nanoribbons show highly nonlinear 
Weber-Fechner like response to static edge potentials. 
These results were obtained from analytical considcra- 
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FIG. 9. (color online) Dependence of magnetization M and 
change in energy AE (referred to the undoped anti-symmetric 
case) as a function of doping S per edge atom. Blue lines cor- 
respond to the symmetric case, while the red lines represent 
the anti-symmetric case. Both M and A are quadratic func- 
tion of the doping 5 for the symmetric case, while they are 
linear for the anti-symmetric case. 



tions of the continuum Dirac equations for graphene. 
These analytical results were validated qualitatively and 
quantitatively with full tight binding calculations. It is 
not often that such a highly nonlinear response of even 
a noninteracting system has been obtained analytically, 
and this is a rare example of a quantum system that 
shows Weber-Fechner like response found in sensory or- 
gans like the eye and ear. It is also interesting to com- 
pare graphene nanoribbons to the well known example 
of a p — n diode (Schockley equation) 2 ^ where the diode 
current depends exponentially on the voltage divided by 
the thermal voltage fcsT/e (/cB~Boltzmann constant, e 
electron charge). In the present case, the Weber-Fechner 
response is an intrinsic zero-temperature property of the 
zigzag edge graphene nanoribbon. 

Our work also suggests many possible uses of graphene 
nanoribbons with applied edge potentials. In the case of 
the symmetric Zeeman potentials, one obtains a system 
with nearly linearly dispersing particles in ID. The veloc- 
ity of these particles can be controlled by tuning the edge 
potential and the width of the nanoribbons. Antisym- 
metric edge potentials provide a route to produce tunable 
ID gapped systems. Furthermore antisymmetric Zeeman 
potentials produce half-metallic cdges i 29 ' 30 Clearly, these 
arc fruitful avenues for further investigation. It will be 
particularly interesting to realize such a physical system 
using cold atom optical lattice o 31 ' 32 where the kinetic en- 
ergy, interactions, and, hopefully soon, temperature can 
be tuned. 

It will be interesting to investigate how physical ef- 
fects not considered in detail in this work affect the non- 
linear response. Primary among them are the second 
neighbor hopping^ 3 - 1', temperature effects, disorder and 
interactions. In particular, the Weber-Fechner like be- 
haviour will require very low temperatures (compared to 
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the energy aeo/W). A simple estimate for the temper- 
ature required to observe Weber-Fechner like response 
can be obtained by comparing the gap scale in (|37|) to 
the temperature. For a ribbon of width ~ 10a, Weber- 
Fechner response should be observable for edge poten- 
tials of order 10 _2 eo — 10~ 4 eo for temperatures below 
~ 10~ 2 eo — 10~ 4 eo. Taking eo ~ 2eV, we estimate this 
temperature to be below 10K. 

The Weber-Fechner like response is likely to be pro- 
tected against weak potential disorder due to the chiral 
nature of the edge states. The most interesting effects 
arise due to the presence of interactions among the elec- 
trons and associated Luttinger liquid effects^ In this 
case, the edge potentials can be "self generated" and can 
lead to interesting magnetic properties of the system. A 
detailed study of these effects has been completed 2 ^ and 



will be reported elsewhere. This study also suggests that 
Weber-Fechner like response is obtained even when the 
edge potential is spread over three/four rows of atoms in 
proximity of the edge. 
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